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^ ' Abstract 

: 

c/3 ' The geometric structure of the Rh^ cation is investigated using a combination of far-infrared 

S . multiple photon dissociation spectroscopy and density functional theory (DFT) calculations. The 

', energetic ordering of the different structural motifs is found to depend sensitively on the choice of 

I pure or hybrid exchange functionals. Comparison of experimental and calculated spectra suggests 

^ ■ the cluster to have a close-packed, bicapped octahedral structure, in contrast to recent predictions 

^ : 

■rj" . of a cubic structure for the neutral cluster. Our findings demonstrate the importance of including 

in ■ 

•n I some exact exchange contributions in the DFT calculations, via hybrid functionals, when applied 

a^ ■ 

• to rhodium clusters, and cast doubt on the application of pure functionals for late transition metal 

o\ : 

. clusters in general. 
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The study of transition metal (TM) clusters offers the opportunity to probe the funda- 
mental physics involved in the transition from atomic to bulk properties |l| and a means 



wide range of cluster properties, includin 
polarizabilityj?] and ionization potential 



to better understand potentially tractable model systems for supported catalysts^]. A 



^, magnetic moments s], |6| , electric 



ig^ reactivity 

9j have been found to depend sensitively, and 
non-monotonically, on cluster size. These size-effects reflect the complex evolution of the 
electronic and geometric structures and yet, the structures of most of these clusters are not 
currently known. 

Although knowledge of cluster structures is vital for developing a deeper understanding 
of the observed properties, experimental structure determination in the gas phase remains 
difficult. Photoelectron spectroscopy has been applied to anionic clusters and has yielded 
detailed information on the size-dependent evolution of electronic structures 10|]. Only in the 



last few years, however, has it become possible to probe more directly the geometric struc- 
tures of transition metal clusters in the gas phase, e.g. by measuring their ion mobilities 
by trapped ion electron diffraction 12|] , or far-infrared multiple photon dissociation spec- 
troscopy (FIR-MPD)Q. 

Theory has been widely used to predict the geometric and electronic structures of clus- 
ters. Electronic structure theory is a vital aid to experimental efforts to determine clus- 
l, [isl, Q] and to help understand the effects of structure on cluster 



ter structure 



properties [15 



12 



161]. The large size of TM clusters, both in terms of the number of electrons 



which must be treated and the complexity of the potential energy surfaces (PES), generally 
restricts practical calculations to the realm of DPT. Application of DFT to mid-sized TM 
clusters requires caution in the interpretation of results; not only because of the approxima- 
tions inherent to all contemporary exchange- correlation functionals, but also because there 
is a lack of benchmark data from high-quality multi-reference electronic structure theory, 
although data does exist for rhodium clusters as large as the pentamer 17|]. Most contempo- 
rary functionals were not designed with TM clusters in mind and it is not clear if one single 
exchange-correlation functional is to be preferred for all TM clusters. 

Rhodium clusters present many of these challenges and opportunities and, as a con- 



This work repre- 



sequence, have been studied experimentally in some detail 
sents the first direct measurements of their structure. Recent DFT calculations have pre- 
dicted a range of unusual putative global minimum and low-energy structures for rhodium 



2 



clusters jisl. 



20| 
14 



2ll| . which are not based on the close-packed motifs identified for 
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many TM clusters 
and a cubic structure for Rh 
clusters 



18 



20 



3 and Rh+Q 
201 ] ■ Cubic motifs have also been reported for larger 



These include a trigonal prism for Rhg 



2l|. Cube-based structures have been rationalized by strong (i-orbital character 



in the bonding, favoring 90° bond angles and eight-center bonding[20]. These structures 
are qualitatively different to the close-packed, polytetrahedral structures reported for most 



cluster sizes by Futschek et al. 



24| . As the computational methods used in these studies 



are rather similar, all using DFT, it is not clear if the differences in the results are due to 
the range of structures and motifs considered for each cluster size or due to differences in 
the details of the calculations, such as the exchange-correlation functionals used. 

Cubic structures have also been predicted for other platinum group metals including 
ruthenium 
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261]. For this metal, Wang and Johnson |2 71] have investigated in detail the 



effects of using pure and hybrid functionals on the favored geometry of small (n < 4) 
clusters. They reported that while pure functionals predicted a square planar isomer as 
global minimum for Ru4, hybrid functionals favored a tetrahedral isomer. They attributed 
the difference to changes in the relative energy of the s- and (i-orbitals in the ruthenium 
atom. The proximity of rhodium and ruthenium in the periodic table suggested similar 
effects may be observed for rhodium clusters, supported by the recent results of Sun et al. 



2l|. Experiments which allow a direct comparison between the experimental results and 



calculated properties therefore provide a vital means to test and benchmark the theoretical 
methods used to study these systems. Here, we focus on structure determination for the 
archetypal cubic cluster, Rhg , by gas-phase vibrational spectroscopy. An eight-atom cluster 
is the smallest cluster which can form a complete cube and therefore provides an ideal 
candidate system to test the competition between cubic and close-packed geometries. IR 
spectroscopy is particularly sensitive to symmetry through its selection rules which, for an 
8-atom cube, predict a single triply degenerate IR-active mode of ti„ symmetry. 

FIR-MPD allows the measurement of size-specific IR spectra of clusters in the gas phase, 
which can be compared to the results of calculations. The technique has been described 
in detail previously [13]. As implemented here, argon-tagged rhodium cluster cations are 
formed by laser ablation of a rotating rhodium metal target. The resulting plasma is cooled 
by a mixed pulse of helium and argon (ca. 0.3%) and carried down a cryogenically-cooled 
clustering channel (173 K) before expansion into vacuum. The resulting cluster beam passes 
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through a skimmer before entering the FIR laser interaction region, where it is intersected by 
a counterpropergating beam from the Free Electron Laser for Infrared experiments (FELIX). 
By recording time-of-flight mass spectra alternately with and without FELIX irradiation, 
the depletion of Ar-tagged clusters can be measured as a function of IR wavelength. From 
this, the size-specific absorption spectra are obtained as described before [13,]. 

The experimental FIR-MPD spectrum of RhgAr"*" is shown in the lower panel of Figure 
[U It appears relatively simple, having an intense peak at 206±1 cm~^ and a smaller peak 
at 250±1 cm^^. The spectrum of RhgAr^ is very similar, showing no large changes in peak 
position or intensity with the degree of argon coverage. This indicates that the Ar acts as 
spectator and is not influencing the cluster structures. This is consistent with the findings, 



e.g. for the larger cobalt clusters 



28|. 



In order to determine the structure of the clusters, the experimental FIR-MPD spec- 
tra are compared with calculated vibrational spectra. A thorough search of the PES has 
been performed in order to identify the important, i.e. low-energy, geometric structures 
of the cluster. We have used a two-stage approach; first employing basin-hopping (BH) 



simulations! 



to locate candidate structures, and second by refining these candidate struc- 



tures. In the first stage, we used BH in two different implementations; the first using the 



Sutton-Chen[30i] mode 
density approximation 



potential, and the second based on the PES described by the local 
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32| . In addition to the structures found in our BH simulations. 



n 



a range of previously reported structures were also investigated 19|, l20|, l24] • In the second 
stage all the candidate structures were further optimized using DFT, performed with the 
Gaussian03 package 331] . Both the pure PBE 3^ and the hybrid PBEl (25% Hartree-Fock 
exchange) 3^, |35| exchange-correlation functionals were used, with the lanl2dz relativistic 
effective core range of spin-multiplicities were considered in 

order to determine the favored spin-state for each isomer and functional. The cluster sym- 
metry was not constrained. All the optimized structures featured Ci symmetry, showing 
small distortions from their higher-symmetry counterparts. Analytic frequency calculations 
were performed to ensure the structures to be true local minima at each level of theory and 
to provide comparison with the FIR-MPD spectra. To aid comparison, the calculated stick 
spectra were broadened by a Gaussian line shape function with a full width at half maximum 
height of 6 cm~^, no scaling of the calculated frequencies was applied. The calculations were 
performed without explicit consideration of the argon tagging atoms as the experimental 
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TABLE I: Relative energies and favored spin-multiplicities of different isomers of Rhg^ at pure and 
hybrid levels of theory. ** Structure collapsed during optimization. 
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spectra were not found to depend significantly on the degree of argon coverage. 

The results of the calculations are summarized in Figure [2] which shows the most im- 
portant geometric structures and in Table [I] where the relative energies of the favored spin 
multiplicity of each isomer at the pure and hybrid levels are listed. We find there to be a 
significant difference in the energy ordering of the isomers at the two levels. At the pure 
level, the cube is the lowest-energy structure followed by the bicapped octahedron, broadly 
in agreement with the results reported by Bae et al. for the neutral clusters [lol. At the 
hybrid level, however, the bicapped octahedron is favored, followed by the bicapped trigonal 
prism. The cube isomer is supported at the hybrid level, but is relatively high in energy 
(0.92 eV). The favored spin-multiplicities for each isomer are generally higher at the hybrid 
level, a point previously noted by Wang and Johnson 27| . This is particularly evident for 
the cube, for which the pure functional favors an octet while the hybrid favors a 12-tet. A 
large number of unpaired electrons is generally in agreement with the experimental finding 
of magnetic moments of up to 0.8±0.2 per atom for small neutral rhodium clustersjsl. 
For ruthenium, significant differences between pure and h ybr id DFT emerge in the predicted 
energetic ordering of the occupied orbitals of the atom 27|. However, for rhodium we do 
not observe such changes, making it more difficult to explain the underlying cause of the 
differences in the pure and hybrid calculations. 

In Figures [H and [3] we compare the experimental FIR-MPD spectrum of RhgAr"*" with 
the calculated spectra from pure and hybrid calculations respectively. Comparison of the 
experimental spectrum with calculated spectra at the pure level show poor agreement with 
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the spectrum of the cubic isomer. For this, three rather closely spaced IR-active modes are 
predicted at ca. 205 cm~^, but no features in the region of 250 cm~^. 

A similar comparison with the spectra calculated at the hybrid level shows the exper- 
imental spectrum to be well reproduced, both in position and relative intensity, by the 
lowest-energy bicapped octahedron isomer. This has an intense feature at 205 cm~^ and a 
weaker feature at 266 cm~^, blue-shifted by 15 cm~^ compared to the experimental spec- 
trum. It is notable that among the different isomers at the pure level the best match to 
experiment is also provided by the bicapped octahedron, though in this case extra features 
are predicted at low wavenumber. The calculated spectra of the other isomers at the hybrid 
level do not provide a good match to experiment, having too few or too many features. The 
spectrum of the cubic isomer is calculated to have a triply degenerate feature of very low 
oscillator strength at 234 cm~^, half way between the two features observed experimentally. 

The combination of observed vibrational spectra and calculations based on the hybrid 
PBEl functional strongly suggest that the RhJ cluster favors a close-packed, bicapped oc- 
tahedron structure and not a cubic structure as has been previously reported, and indicated 
by our own DFT calculations based on the pure PBE functional. In the absence of high-level 
multi-reference benchmark calculations, our findings cast doubt over the suitability of pure 
functional when applied to small/mid-sized late TM clusters. Our evidence, while specific 
to Rh^, suggests in general that the expectation of open structures such as those based on 
cubic motifs in a range of late TM cluster sizes may be misconceived. Further evidence, 
comprising close comparison of experiment and theory for other sizes of Rh clusters, will be 
required to resolve this question. 
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FIG. 1: Predicted IR spectra of Rh^ from pure (PBE) DFT calculations compared to the ex- 
perimental FIR-MPD spectrum of RhgAr^, a is the experimental IR cross section. The labels 
correspond to the structures shown in Figure [2j 
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FIG. 2: Geometric structures of low-energy isomers of Rh^ identified in DFT calculations: a) 
bicapped octahedron (bc-oh), b) bicapped trigonal prism (bc-tp), c) cube, d) diamond prism 
(diam), e) square antiprism (sq-ap). 
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FIG. 3: Predicted IR spectra of Rh^ from hybrid (PBEl) DFT calculations compared to the 
experimental FIR-MPD spectrum of RhgAr^. The labels correspond to the structures shown in 
Figure [3 
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